A new method for analyzing second-order phase transitions applied to the 
ferromagnetic transition of a polaronic system 
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A new method for analyzing second-order phase transitions is presented and applied to the pola- 
ronic system Lao.7Cao.3Mn03. It utilizes heat capacity and thermal expansion data simultaneously 
to correctly predict the critical temperature's pressure dependence. Analysis of the critical phenom- 
ena reveals second-order behavior and an unusually large heat capacity exponent. 

PACS numbers: 75.40.Cx, 65.20.+W, 75.47.Lx, 75.30.Kz 



Investigations of phase transitions are important for 
our fundamental understanding of condensed matter sys- 
tems and of wider interest because these systems can 
serve as environments for the study of topological defects 
Q . The phase transition in manganese oxides exhibiting 
colossal magnetoresistance (CMR) is particularly fasci- 
nating since competition among the charge, lattice, and 
spin degrees of freedom leads to the CMR effect 0. 
Double-exchange interactions among the magnetic ions 
and electron-phonon coupling via the Jahn- Teller distor- 
tion play essential roles, especially near the ferromagnetic 
(FM) to paramagnetic (PM) phase transition where the 
metal-insulator transition occurs and the CMR effect is 
the largest. Defects known as magnetic polarons form 
well above the critical temperature T c and increase in 
density as the temperature is lowered through T c 0, [j| . 
The complexity of this phase transition has led to uncer- 
tainty as to its thermodynamic characterization 0. 

The first-order phase transition line, along which two 
thermodynamic phases coexist on a phase diagram, is 
described by the Clausius-Clapeyron equation. Across 
this line the first derivative of the free energy is discon- 
tinuous. This results in the divergence at T c of quanti- 
ties, such as the molar heat capacity at constant pres- 
sure Cp or the thermal expansion coefficient fi, but no 
sign of divergence occurs prior to approaching T c . In 
general, this generic feature is not compromised by fi- 
nite system size because the correlation length remains 
finite at T c and is typically much smaller than the sam- 
ple size. In contrast, the second-order (continuous) phase 
transition occurs with diverging correlation length as T c 
is approached, indicating a scale-free phenomena. It is 
possible to generalize from first-order transitions with 
discontinuous first derivatives of the free energy to nth- 
order transitions with discontinuous nth-derivatives 
Indeed, the Ehrenfest equation 
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is a generalization of the Clausius-Clapeyron equation to 
second-order transitions (v is the molar volume, ACp 



and AO = 3A/i are the jumps in Cp and the volume 
thermal expansion coefficient fi at T c , and P is the pres- 
sure). Interestingly, the only phase transition known to 
exhibit distinct jumps in Cp and ft and obey Eq. (1) is 
the normal-superconductor transition . For most sys- 
tems undergoing continuous phase transitions, instead of 
jumps in Cp and £7, critical behavior is observed. More 
specifically, near T c , Cp (with t = (T — T c )/T c ) is domi- 
nated by Cp ~ |t| _Q± . The exponents a± are called the 
critical exponents above (+) and below (— ) T c . Appar- 
ently, if Cp or Q exhibit critical behavior, the quantities 
ACp and Af2 in Eq. (1) are not well defined. One result 
of this letter is a substitute for Eq. (1), that is widely 
applicable for calculating dT c /dP using C p and Q data. 

With T and P as independent variables, the thermo- 
dynamic potential per mole $ = u — TS + Pv is used to 
derive the needed thermodynamic equality (u and S are 
the internal energy and entropy per mole). Let T C (P) be 
the critical temperature at pressure P; it marks a phase 
transition line on the (T, P) plane. We assume, within 
the pressure range of interest, that T C (P) can be inverted 
into P{T C ). The molar entropy at T C (P), S(T C ,P(T C )), 
abbreviated by S(T C ), marks a transition line on the 
(T, S) plane as well. Given a small temperature shift 
8 away from a fixed T c and starting at (T c + 5, P(T C )) on 
the (T, P) plane, we move an infinitesimal amount par- 
allel to the transition line P(T C ). It is easy to show that 
such movement also moves the point (T c + 8, S{T C )) on 
the (T, S) plane parallel to the line S(T C ). The displace- 
ments dS, dT, and dP are related through 



TdS = CpdT - vTVldP. 



This immediately leads to 
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which holds for T = T c + 8 above and below T c ; (dS/dT) c 
and (8P/dT) c are slopes of the phase transition lines at 
fixed T c on the (T, S) and (T, P) planes, respectively. 
Eq. (3) implies that, near T c , one can superimpose C* = 
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C p — a — bT with 7] Til after rescaling of Tfl; v is treated 
as a constant since it changes only slightly, < 0.1%, near 
T c . This relation also indicates that if Cp diverges, Tfl 
diverges with the same exponent. 

We investigate a CMR oxide for which it has been dif- 
ficult to determine whether the phase transition is first 
or second order 0, S H El El E3- Polycrystalline 
Lao.7Cao.3Mn03 was prepared by mixing stoichiomet- 
ric amounts of La(N0 3 ) 3 -6H 2 0, Ca(N0 3 ) 2 -4H 2 0, and 
C4HgMn04-4H20 in distilled water and a 50 mol% excess 
of citric acid and ethylene glycol. The solution was heated 
at 120°C, stirred until a gel formed and then dried. The 
organic material was oxidized 24 h at 500° C. The powder 
was ground with an agate mortar for 30 min, heat treated 
30 h at 1000°C, ground for 30 min, and heat treated 30 
h at 1100°C. Finally, the powder was ground for 30 min, 
pressed into pellets and reacted 30 h at 1200°C. The den- 
sity was 5.47(2) gm/cm 3 (90% of theoretical density) and 
x-ray powder diffraction confirmed the single-phase na- 
ture. Four-probe electrical resistivity and ac susceptibil- 
ity (500 Hz) were measured under hydrostatic pressure 
in a Fluorinert medium; a manganin sensor determined 
pressure at 295 K. A fused quartz capacitive dilatometer 
was used to observe the linear thermal expansion with a 
sensitivity in Al of 0.1 A; for the 2.6 mm long sample, 
the relative sensitivity is about four orders of magnitude 
better than diffraction methods. Heat capacity and dc 
magnetization were measured with a Quantum Design 
PPMS. All data were obtained on the same sample. 

C P and S {S{T)= J (C P /T)dT) are shown in FigH^a). 
The anomaly in Cp is associated with the PM to FM 
transition. To obtain S associated with the magnetic 
transition, Smag, a polynomial fit from 80 K to 325 K 
(excluding the region 200 K < T < 270 K) was sub- 
tracted; Smag (upper inset of Fig. ^a)) reveals a 
smooth change at T c . The entropy value 3.2 J/mol.K at 
T = 300 K, agrees with prior reports 

[HEi. Th e ther- 
mal expansion coefficient n(T) is illustrated in Fig. db); 
H was determined by taking a point-by-point derivative of 
AI/Iq (shown near T c in the inset), it reveals an anomaly 
at T c as well. Straight lines are drawn through Cp and 
/i above and below T c illustrating the form of an ideal 
second-order phase transition jjj. This provides AC p — 
-7.45(5) J/mol.K and Att = 3A^ = -9.60(12) x 10" 6 
K _1 . Eq. (1) is applied to estimate the pressure deriva- 
tive dT c /dP = 11.5(2) K/GPa using T c = 249.5 K and 
v = 3.57xl0~ 5 m 3 /mol. Magnetic susceptibility at 2000 
Oe, n(T) and C P {T) all reveal T c = 249.5 K. 

To apply Eq. (3), the overlap between C* Q and A/iT 
is maximized with a = 59.5 J/mol.K, b — 0.128 J/mol.K 2 
and A = 7000 ± 300 J/mol.K. The proportionality be- 
tween C* and fiT (Fig. 2) suggests: (1) S must be con- 
tinuous and differentiable along (and in the vicinity of) 
the phase transition line and (2) v at a given pressure 
is continuous across T c ; therefore, the phase transition 
is continuous (second order). Furthermore, the quantity 
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FIG. 1: (a) Molar heat capacity Cp, molar entropy S, and 
magnetic entropy Smag (inset) versus T. (b) Expansion co- 
efficient /j, and Al/lo (inset); lo is the length at 300 K. Lines 
denote the jumps for an ideal second-order phase transition. 

dT c /dP is now given by 

§-(£);' -^- im (7)K/Gpa - <4 » 

In Fig. 2, Cp* displays more rounding of the peak than 
\x\ this is associated with the heat pulse which changes 
the sample temperature by AT (~ 2.5 K near T c ) at each 
measurement temperature. This leads to an averaging 
effect ^3 that is absent in the thermal expansion data 
which is acquired closer to equilibrium (warming rates < 
8 K/h). In addition, there is a smaller effective sample 
size in the Cp data because the heat pulse cannot warm 
the sample uniformly and simultaneously. 

Before comparing dT c /dP in Eq. (4) with experiment, 
critical exponents are determined. This usually requires 
a background subtraction to enlarge the range of validity 
of power-law behavior. Although more complicated back- 
ground subtractions might, in general, be necessary, we 
allow only two constant shifts A± in the expression of /xT. 
That is, we assume that when T > T c (T < T c ), \iF can 
be fitted by the expression A + (_) + (i3 + (_) /a) \t\~ a . Mak- 
ing a grid for -12 < A± < 12 and 248K < T c < 251K 
and plotting \t\ versus 7000/iT — A± on a log- log scale, we 
choose the set of A± and T c that maximize the power-law 
fittable temperature range and minimize the deviation 
within the fitted range (220 K < T < 270 K). Using the 
optimal A± and T c , Cp - A± and 7000//T- A± versus t 
are plotted in Fig. 3. The ability to obtain a power-law 
fit above and below T c with the same exponent agrees 
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FIG. 2: Molar heat capacity after subtracting a linear term, 
Cp* (open symbols), and 7000//T (closed symbols) versus T 
illustrating the scaling suggested by Eq. (3). 




FIG. 3: Cp - A± (open symbols) and 7000^T - A± (closed 
symbols) versus t = (T — T c )/T c on a log-log scale. 

well with the one-parameter scaling theory of continu- 
ous phase transitions. The smaller reduced temperature 
range for T > T c indicates that the fluctuation-induced 
ordered domains grow to macroscopic size only when T 
is close to T c . This implies that the prefactor for the cor- 
relation length £ ~ g±\t\~ 1 ' is smaller on the high tem- 
perature side (g+ < gJ). The exponent a± — 0.93(8) 
is obtained whose magnitude is significant since a± > 1 
would lead to a divergence in the integrals of Cp and /i 
near T c 0. This explains why this continuous (second- 
order) phase transition, evidenced by the superposition 
of C* and 7000/xT, has often been characterized as first 
order or nearly first order. The same measurements and 
analysis on a second sample of La .7Cao.3Mn03 yielded 
identical results for a± |l8|. 



FIG. 4: Electrical resistivity and ac susceptibility (upper 
inset) versus T under pressure P. Lines illustrate T c determi- 
nation. Lower inset shows T c and Tmi versus P. 

To test Eqs. (1) and (4), dT c /dP was measured. Fig. 
0] shows the electrical resistivity p(T,P). Pressure de- 
creases p and shifts the metal-insulator transition tem- 
perature Tmi upward. In the left inset of Fig. 0] the 
pick-up coil signal is shown and lines illustrate the man- 
ner in which T c was determined. The (T, P) phase di- 
agram is shown in the lower inset. These results re- 
veal that dT MI /dP = dT c /dP = 16.5(3) K/GPa [3, 
which is 43% larger than the value of 11.5(2) K/GPa 
calculated by Eq. (1), but in better agreement with 
15.3(7) K/GPa calculated by Eq. (4). The wide range 
of dT c /dP values (12 to 22 K/GPa) 0, for this com- 
position (x = 0.30 and 0.33 values are noted here) il- 
lustrates that all quantities appearing in Eqs. (1) and 
(4) must be determined on the same specimen in order 
to conduct a meaningful analysis. 

As noted above, agreement has not been reached 
regarding how to characterize the phase transition of 
Lao.7Cao.3Mn03 (and compositions close to this). Heat 
capacity measurements on Lao.67Cao.33Mn03 were un- 
able to determine the critical exponent indicating that 
the phase transition was first order; however, critical ex- 
ponents for Lao.6Cao.4Mn03 were obtained, which led 
to the assertion that it defines a tricritical point in the 
phase diagra m 1 121 . Our results together with those 
of Kim et al.|l2j, indicate that further investigations 
are needed to determine the compositional dependence 
of a±. An analysis (La .65Ca .35MnO 3 ) using the 
Clausius-Claperyon equation suggested that the phase 
transition was first order. We have also applied this 
method and found dT c /dP in agreement with our mea- 
surements; however, this is fortuitous since the Clausius- 
Claperyon equation is valid for transitions where S and 
v make discontinuous changes at T c , which is clearly not 
the case (see insets of Fig. 1). In comparing our study 
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to others, it is essential to consider: (1) the scaling be- 
tween C* and fiT, (2) the use of both in obtaining the 
critical exponent and (3) a± is close to 1 resulting in 
near— divergences of Cp and fj, which clarifies why the 
transition has often been identified as first order. 

Establishing the continuous (second-order) nature of 
the phase transition has important implications. Above 
T c , magnetic defects (polarons) form which grow in den- 
sity as the sample is cooled through T c Our re- 
sults suggest a divergent magnetic correlation length at 
T c only if the magnetization is the correct order parame- 
ter in CMR systems. Recent neutron scattering measure- 
ments suggest a large, but probably finite magnetic 
correlation length which was argued to suggest first-order 
behavior. However, since our thermodynamic analysis 
does not require direct measurement of the correlation 
length, it stands correct without prior specification of 
the order parameter. 

The extremely large value of a± — 0.93(8) implies 
smaller exponents for the magnetization and magnetic 
susceptibility if the exponent identity a + 2/3 + 7 = 2 
holds. The smaller /3 and 7 exponents would indicate 
weakening of the effective magnetic coupling as T c is ap- 
proached and are likely a direct consequence of compe- 
tition among charge-lattice-spin coupling. This suggests 
that the magnetization, alone, may not be the appropri- 
ate order parameter. However, as a cautionary point, it 
is noted that the thermodynamic relations impose only 
the inequality a+2/3+7>2asa restriction |20j. It is im- 
portant to consider that a± = 0.48(6) and 0.05(7) were 
determined for Lan fiCap 4Mn0.s and Lao.7sSro.25Mn03, 
respectively [l2j, |21| . Their smaller ct± values (or larger /3 
and 7 values) can be attributed to weaker charge-lattice- 
spin coupling as ascertained by smaller CMR effects rel- 
ative to La .7Cao.3Mn03 jl^. 

Finally, it is interesting to consider topological defect 
experiments. In superfluid helium, pressure is applied 
to move out of the superfluid state, then quickly released 
while measuring the defect concentration (superfluid vor- 
tex density) as a function of time |2flj . Similar experi- 
ments in CMR systems could be envisioned where pres- 
sure rapidly moves the system into the FM state while 
the concentration of magnetic defects is observed versus 
time; such experiments are simplified by a large dT c / dP, 
strong magnetic signal, and high temperatures. 

In summary, a new method for investigating continu- 
ous phase transitions that uses heat capacity and thermal 
expansion data was presented. Its application to the FM 
phase transition of the CMR oxide Lan.7Can.3Mn03 has 
revealed an extremely large heat capacity exponent. 
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